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We compute the finite temperature induced fermion number for fermions coupled to a static 
nonlinear sigma model background in (2 + 1) dimensions, in the derivative expansion limit. While 
the zero temperature induced fermion number is well known to be topological (it is the winding 
number of the background) , at finite temperature there is a temperature dependent correction that 
is nontopological - this finite T correction is sensitive to the detailed shape of the background. At 
low temperature we resum the derivative expansion to all orders, and we consider explicit forms of 
the background as a CP instanton or as a baby skyrmion. 
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The phenomenon of induced fermion number arises due to the interaction of fermions with nontrivial topological 
I backgrounds (e.g., solitons, vortices, monopoles, skyrmions), and has many applications ranging from polymer physics 
to particle physics jl], g, [|, |], 0, || . At zero temperature, the induced fermion number is a topological quantity 
(modulo spectral flow effects), and is related to the spectral asymmetry of the relevant Dirac operator, which counts the 
difference between the number of positive and negative energy states in the fermion spectrum. Mathematical results, 
such as index theorems and Levinson's theorem, imply that the fractional part of the zero temperature induced 
, fermion number is topological in the sense that it is determined by the asymptotic properties of the background fields 
£NJ ' H [To], [n], [l3|, . This topological character of the induced fermion number is a key feature of its application in 
certain model field theories jl^, |lG. 17, At finite temperature, the situation is very different - the induced fermion 



number is generically nontopological, and is not a sharp observable 19 |2(i| , plf. Several explicit examples of kinks 



^2], 23, |2J and sigma models p0| in (1 + 1) dimensions, and monopoles ]19[ |26| , p7| , [28| in (3 + 1) dimensions, 
have been analyzed in detail. In this paper we compute the finite temperature induced fermion number for fermions 
| coupled to a static nonlinear sigma model background in (2 + 1) dimensions. 

. At zero temperature, the induced fermion number for fermions coupled to a static nonlinear sigma model background 
in (2 + 1) dimensions has been studied extensively |2^, [30| |3l], pS2fl . The T = induced fermion number is equal to 
the winding number of the sigma model background field, and this result may be interpreted in terms of a topological 
current density. This system has applications both in condensed matter physics |n], where the sigma model 
provides a phenomenological model for a two dimensional Heisenberg ferromagnet, and in particle and nuclear physics, 
where the sigma model can be used as a "baby Skyrmion" model P9| , |30| , |33| , |35f that mimics many properties of the 
. £h \ (3 + 1) dimensional Skyrme model for baryons p6| . One motivation for this paper is to compute the induced fermion 
number at nozero temperature and to understand in detail the origin of the nontopological T dependent contributions, 
$H ' as similar effects will occur in the (3 + 1) dimensional Skyrme case. 

In Section II we define what is meant by finite temperature induced fermion number, and indicate how it can be 
computed. In Section III we use the derivative expansion to compute the finite temperature induced fermion number 
for a nonlinear er-model background in (2 + 1) dimensions. At low temperature we resum the derivative expansion to 
all orders, and we consider explicit forms of the background as a CV 1 instanton or as a baby skyrmion. We conclude in 
Section IV with some comments concerning the relation of our results to the well known T = results, and concerning 
the possible extension to (3 + 1) dimensional Skyrme models. 



II. FINITE TEMPERATURE INDUCED FERMION NUMBER 

The induced fermion number is an expectation value of the second quantized fermion operator N = ^ J dx [\&' , 
For a given static classical background field configuration, the fermion field operator ^ can be expanded in a complete 
set of eigenstates of the Dirac Hamiltonian H. Expectation values of N can then be computed. At zero temperature, 
the fermion number is a vacuum expectation value (N)q = (0 |iV| 0) , and is related to the spectral asymmetry of the 
Dirac Hamiltonian 

(N)q — (spectral asymmetry) 
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dEa(E) sign(£) (1) 



Here a(E) is the spectral function of the Dirac Hamiltonian H : 



= — Im Tr ( — -1—: - 1 (2) 



7T \H-E-ie 

At nonzero temperature, T, the induced fermion number is a thermal expectation value : 

Tr (e~^ H N) 



Tr (e-^) 

-1 r dEa(E) tanh(^) (3) 



where /? = y . Notice that this finite temperature expression (^|) reduces smoothly to the zero temperature expression 
(1) as (3 — > oo. In fact, the nonzero temperature provides a physically meaningful and mathematically elagant 
regularization of the spectral asymmetry. 

The spectrum of the fermions is independent of the temperature in our approximation of a static classical background 
field. All information about the fermion spectrum, in the presence of the background, is encoded in the spectral 
function a(E) defined in @. Thus, knowledge about the spectral function a(E) is the key to evaluating ([!]) or (||). 
Actually, to compute the induced fermion number (either at T — 0, or T > 0) one only needs the odd part of the 
spectral function, as is clear from ([[]) and (||). Thus, the calculational problem is to find the odd part of <j{E). 

A more physical interpretation of the finite T fermion number (N)t is provided by noting that (N)t in (|^) 

separates naturally into a T — piece and a finite T correction 



(N) T = {N) + dEa(E)sign(E)n(\E\) , (4) 

J —OO 

where 

n{E) = ^h~i (5) 

is the Fermi-Dirac distribution function, and we have used the simple fact that tanh(^) = 1 — 2n{E). 

The first term, (iV)o, in (^) is known to be topological in the sense that it is invariant under small local deformations 
of the background which do not change the asymptotic boundary values of the background. Physically, this topological 
term corresponds to the vacuum polarization response of the system to the presence of the background. The second 
term in (Q), which is the finite T correction term, is generically nontopological in the sense that it is sensitive to 
local deformations in the background. Physically, this term corresponds to the plasma response of the system to the 
presence of the background. It incorporates effects beyond the vacuum response, and so is more sensitive to the entire 
single-particle fermion spectrum. The caveat "generically" has been included here because for certain very special 
backgrounds the finite T correction is itself topological. This happens when the Dirac Hamiltonian has a quantum 
mechanical supersymmetry relating the positive and negative parts of the spectrum pp| , [2l| ] . Several explicit examples 
(in (1 + 1) and (3 + 1) dimensions) of this distinction between topological and nontopological contributions have been 
studied in detail already [|l9[ An example in (2 + 1) dimensions where the finite T correction is topological is 

the case of a background flux string |2(| |2^| . In this paper we study a (2 + 1) dimensional example in which the 
finite T correction is nontopological. 

III. THE NON LINEAR SIGMA MODEL IN (2+1) DIMENSIONS 

Consider fermions coupled to a static er-model background in (2 + 1) dimensions. 

C = ^{i^d^ — rah ■ r)V> (6) 

The fermions are in the defining representation of SU(2), and the static cr-model background is represented by a 
3- vector (in internal space) h(x), which is constrained to take values on S 2 : 

h\x) = 1 (7) 
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In the Lagrangian (||), the r° are Pauli matrices, and to is a mass parameter that sets an important energy scale 
of the single-particle fermion spectrum. Later in the paper we will consider specific profiles for n(x). At T = the 
quantization of fermions in this background is well known, and leads to an induced fermion number |3(], ||lj 

(N) = -!- [ d 2 x e abc e l] n a d t n b d,n c (8) 

»7T J 

This induced fermion number, (N)o, is just the integer winding number of the map n : Sft 2 — > S 2 . With the boundary 
condition that h(x) has an angle independent limit at spatial infinity, we can compactify the spatial 5ft 2 to S 2 , in 
which case {N)o gives the integer winding number for maps n : S 2 — > S 2 . Thus, (N)o is topological. 

Invoking Lorentz invariance, one interprets the result (0) as the spatial integral, J d 2 xJ° , of the zeroth component 
J° of a "topological" current density 

J" = — e"" A e abc n a d v h h d x n c (9) 
8n 

This current is manifestly conserved (d^J^ — 0) without use of the equations of motion, and hence is called "topologi- 
cal". In fact, J M may (and generally does) have higher derivative corrections beyond the term shown in ([)]). However, 
these are all total derivatives and hence do not contribute to an integrated quantity such as the zero temperature 
fermion number (N) in (§) |o|. 

At nonzero T, several aspects of this familiar story change. First, as we show explicitly below for this (2 + 1) 
dimensional cr-model case, there are higher derivative corrections to J° that are not total derivatives, and these 
contribute to the finite T correction to the T — induced fermion number (||). Second, at finite T one can no 
longer invoke Lorentz invariance (the thermal bath determines a preferred frame), and so one cannot directly infer a 
topological current density as was done in (||). 

A convenient approach to this calculation is to express (|§|) as a contour integral (using the ie prescription in the 
spectral function) in terms of the Dirac resolvent 

R(z) = Tr(-^-L-) (10) 



Thus we have 



, , 1 f dz _ , 1 v , 0z , , 

(N t = -r / ^-t Tr tanh V 11 

2 r 27TI H — z I 



where C is the contour (— oo + ie, +oo + ie) and (+oo — ie, — oo — ie) in the complex energy plane. For a hermitean 
Hamiltonian H, the poles and branch cuts of the resolvent lie on the real axis, so one has two choices for evaluating 
the contour integral in (|ll]). First, one can deform the contour C around the poles and cuts of the resolvent R(z). 
This approach leads to an integral representation for (N)t- Alternatively, one could deform the contour C around 
the Matsubara poles, z n — (2n + 1) iirT, of the tanh(^) function, which lie on the imaginary axis. This leads to a 
summation representation for (N)t- These two representations are, of course, equivalent. The integral representation 
is just the familiar Sommerfeld- Watson transform of the summation representation. 

Thus, in order to evaluate the finite temperature induced fermion number, (N)t, one needs the odd part of the 
spectral function cr(E); or, equivalently, the even part of the resolvent R(z). No exact expression is known for 
Rcvcn{z) for the interaction described by (|^). Therefore, one needs an approximate technique. Here we use the 
derivative expansion [^6| , in which one assumes that the natural length scale A associated with the background field, 
n(x), is large compared to the Compton wavelength ^ of the fermion field: Am 3> 1. That is, we assume that 

|Vn(x)|<m (12) 

By symmetry, this means that we assume that all derivatives of n are small compared to the fermion mass scale m. 
Since the background is static, these derivatives are all spatial derivatives. The derivative expansion of the resolvent 
(and hence of the induced fermion number) can then be formalized as follows. The Dirac Hamiltonian corresponding 
to the Lagrangian (||) is 

H = -id ■ V + to h ■ t 7° (13) 

Here, (a = 7 7), and for definiteness we choose gamma matrices as follows: 7 = a 3 , 7 = —ia. The even part (in z) 
of the resolvent (|l(]) is 



Tr 



H-z 



I ( " Tr 7 " ~ ) < l4 > 
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This can be systematically expanded in terms of derivatives of the background field as follows. First, write the square 
of the Hamiltonian (|l3|) as 



H 2 



— V 2 + m 2 + im^djfi ■ f 



H$ + V 



(15) 



where H 2 = — V 2 + to 2 is the square of the free Hamiltonian, and the "interaction" is V — irrvy 3 djn ■ f. Next, since 
V involves a factor of m and a derivative of n, by ( |l2| ) it is small relative to the natural scale, m 2 , of H 2 . Therefore, 
the even part of the resolvent in (114) can be systematically expanded in terms of V : 



Tr 



H 



= Tr(AFAFAJ) - Tr(AVAFAFAlf) + Tr(AVAFAFAFAJ) 



(16) 



where we have defined the kinetic part of H as K = — id ■ V, the interaction part of H as I = mh ■ T7 , and the free 
propagator as 



A = 



1 



(17) 



In obtaining ( [l6| ) we have used the simple facts that Tr(r a ) = and Tr(7 n . . .j %an+lr y D ) = 0. The trace in ( fl6| ) 
involves a matrix trace over the Dirac matrices, 7 M = (a 3 ,— it?), and the isospin matrices, r, as well as a functional 
trace over the propagators A. These traces can be done in momentum space or in position space. In position space, 
as we go to higher derivative orders we make use of the following operator identity: 



1 



-V(x) = V{x)- 



1 



-V 2 



2 W(f) • V 



1 



(-V 2 + a 2 ) 2 



(vV(f)) 



1 



+4 (didjV{$))didj- 



1 



-V 2 



-,2^3 



V (V 2 V(x) 



(-V 2 + a 2 ) 3 



(18) 



where the parentheses around the derivatives on the RHS of ( |T8| ) indicate that the gradient operator acts on V(x) 
only. At any given order of the derivative expansion we collect together the required number of derivatives of h{x) , 
as given by repeated application of the above identity, noting that V(x) itself involves one derivative of h(x). This, 
however, becomes rather clumsy for higher order terms in the derivative expansion and in practice it is easier to do 
the computation in momentum space. The leading term in the derivative expansion involves two derivatives, which 
come directly from the two factors of V in the first term of (Eq) : 



Tr 



H-z 



(2) 



= Tr (AVAVAI) 

— TO 

2tt(to 2 - z 2 y 



(2) 



d 2 xe l] e abc h a dih b djh c 



(19) 



To compute the induced fermion number to this (leading) order of the derivative expansion, we evaluate the z integral 
in (11) as a sum over Matsubara modes. This leads to the following contribution to the induced fermion number: 



<JV#> 



E 

l — — DC 

OO 

E 



to T 



27r[m 2 + ((2n+l)7rT) 2 ] 2 
4m 3 T \ 



d 2 xe ij e ahc h a dih h d,n c 



(N)o 



\n— — oo 



[m 2 + ((2n + 1)^T) 2 ] 2 / 

The prefactor in (|20|) is temperature dependent, but in the T — * limit it reduces smoothly to one: 



E 



4m 3 T 



= 1 - 



n— — oo 



2m 



[m 2 + ((2n+l)^T) 2 ] 2 
Here we have used the general low temperature expansion: 

1 m 1 ~'* p T(j3-\) (2mT) 1 -P _ m 



* E 



(((2n + l)7rT) 2 + m 2 )P 2^ T(p) 



mT(p) 



(20) 



(21) 



(22) 
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Thus, in the T — > limit, the leading derivative expansion term ( pp| ) produces the entire zero temperature answer, 
(N) . And at nonzero temperature, this order of the derivative expansion gives the result (^0|), which is topological 
because it is simply the T = result multiplied by a smooth function of temperature. 

To see the nontopological contribution to (N)t we need to go to the next order in the derivative expansion, which 
involves four derivatives of n(x). Recalling that V = im^ J djh{x) ■ f already includes one derivative of h(x), we see 
that at fourth order we can get contributions from the first three terms in (|16[): 



Ti- 



ff 



(4) 



Tr(AV/AVA/) (4) - ^T{AVAVAVAK) (/^, + Tv{AV AV AV AV AI) 



(4) 



\(4) 



(23) 



After a straightforward expansion one finds 



Tr 



1 



H-z 



(4) 



■5z" 



12n (m 2 - z 2 ) 4 



12tt(to 2 - z 2 ) 3 



d 2 x (e ijk e abc n a dih b d ] n c ){d k h d d k h d ) 



d 2 s 



2V 2 (e ij ' e abc h a dih b djh c ) + 3d l (e ij e abc (V 2 n Q ) djft b ri c )] (24) 



There are several important comments to be made about this fourth order contribution. First, the second term 
contains total spatial derivatives, which vanish after doing the x integral. Second, the first term is nontopological 
and cannot be expressed as a total derivative. This term gives a nonzero contribution to the spatial integral, and 
furthermore this contribution (unlike the winding number integral in (|l9|)) depends explicitly on the length scale of 
the background field n(x). Third, to see how this is compatible with the fact that the zero T fermion number is 
topological, we observe that after doing the z integral in (pd|), we obtain 



(4) 



m 2 - 5((2n+ 1)ttT) 2 
12tt ^ \m 2 + ((2n + l)7rT) 2 )l 4 

11— — OO 



E 



d 2 x (e tj t abc h a dih b djn c )(d k n d d k h d ) 



(25) 



The prefactor is T dependent, as was the prefactor in the 2-derivative contribution in (|20|). However, using (|22j), as 
T — > 0, the prefactor of the fourth order term (Eq) behaves as 



m 3 T m 2 - 5((2n + l)7rr ) 2 

12tt f-^ [m 2 + ((2n + 1)ttT) 



214 



rn 
12^ 



3 r 



( 15 



1 



\16m 5 8m 2 T 3 
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16m 5 8m 3 T 2 



/ m 
V967rT 3 



T 

1 - 5 — 

m 



(26) 



so that it vanishes (exponentially) as T — > 0. Thus, at T = the nontopological contribution to the fermion number 
vanishes, at this order in the derivative expansion. Note that in ( fj"l|) the leading constant term 1 survives the T — > 
limit, while in ( p6] ) the leading constant terms cancel, leaving a function that vanishes as T — ► 0. 

Already at this next-to-leading order of the derivative expansion, we have established in ( ]25| ) that the induced 
fermion number acquires a nontopological temperature dependent contribution at finite T. We now consider higher 
orders in the derivative expansion. The next term beyond (|24|) has 6 derivatives and receives contributions from 



each of the first five terms in the expansion (16). While it is systematic to write down all these 6-derivative terms, it 



becomes a lengthy expression due to the many ways of contracting the spacetime and SU(2) indices. This proliferation 
of terms rapidly becomes worse as we go to higher and higher orders in the derivative expansion. To proceed, we 
consider instead the low temperature limit. In this case we find that we can evaluate the leading contribution to 
all orders of the derivative expansion, and rcsum this leading low temperature correction to the zero temperature 
answer (||). This follows the procedure used in ||l9| f or a 1 + 1 dimensional sigma model, and is based on simple 
dimensional analysis. From the general expansion (|16|) we see that terms either involve: (a) an even number of V's 
and one interaction term J, or (b) an odd number of Vs and one kinetic term K. Consider a term of form (a) with 
v vertex insertions of V. Since V has already a derivative of n(x), the total number of derivatives, d, of h(x) in this 
term satisfies 2 < v < d. This term also has v + 1 propagators A. Thus, by dimensional reasoning, this term must 
behave as 



T>((AV)"AJ) 



(<*) ~ A 

(m 2 - z l)(v+d)/2^d 



nix) 



(27) 



where A v d [n(x)] is some functional of h(x) with d derivatives and (y + 1) factors of n(x). On the other hand, a term 
of the form (b) behaves as 



Tr((AV)"AK) 



(d) 



( m 2 _ z 2yu+d-l)/: 



(28) 
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We now consider the T dependence of the prefactors of these types of contributions. First, for d > 4 the constant 
terms in the low T (T -C m) expansion must all cancel at any given order d of the derivative expansion, since the 
final answer must be just the 2-derivative topological term (||) at T = 0. This is a very stringent test of the derivative 
expansion at finite temperature. We have explicitly checked that this is indeed satisfied for d = 4 and d = 6. The 
leading corrections all have the same exponential factor, with a prefactor coefficient that depends on the power of 
the propagator. From (|2^) we see that the dominant prefactor at low temperature occurs when p (the power to 
which the propagator is raised) is as large as possible. This means that for the terms in ( p^ ) and ( p8| ) we must take 
v = d. This means that each derivative just comes from one of the insertions of V. This fact dramatically simplifies 
the derivative expansion (in this low T limit), because the propagator factors can simply be moved out of the trace 
without generating further derivatives (since these would be subleading). We also learn that the terms of the form ( p7| ) 
dominate over those of the form (|2^). These simplifications permit us to go to any order of the derivative expansion. 

To compute the actual form of the derivtive expansion contribution in this low temperature limit, we have to 
calculate ^ +1 [n(x)] in (|27|), which involves doing the (matrix) traces over the Dirac and isospin matrices: 



tr [lV d ] = tr 



m a 



<g> (n(x) ■ f) (m 2 d l h a d 3 n b a 1 a ] 



"Y 



(29) 



where (er 3 , —ia) and f are the Pauli matrices in the Dirac and Isospin spaces respectively. Using the fact that i and 
j can only take values 1 or 2, and that a l a^ — <5 y + i e v a k , the above trace becomes 



m d+1 tr 



{h(x) ■ t) (dih a dih a - e lj e abc dih a djh h a z <g> t c ) 1 



(30) 



Using the binomial expansion for expanding the above power and noting that: (i) Tr (cr 3 ) m+1 = 2 when m is odd and 
zero otherwise, and (ii) (r • a)(r • b) = a ■ b + i(a x b) ■ f , the expression (|3C|) becomes 



4m 



d+ 



1 g(g) 

\m\ 



d/2 

E 



2 ^{d l n a d i n a )^~ m) \v{x)\ 7 
m 1 



m— 1,3,5 

where we have defined the following combinations 

v{x) = e lj dih x djfi , q(x) = e ij e abc n a djfi c 
Finally, performing the functional trace over the propagators we find 

r q (x) 



(31) 



(32) 



Tr((AU) t/ A/)( d) 



27rd(m 2 - z 2 ) 



d z x 



v{x)\ 



{d t h a d t h a + \v(x)\) d/2 ~ {d l n a d % h a - \v{x)\f 2 (33) 



So, in the low T limit (T <C m), the derivative expansion becomes very compact. The leading low temperature 
contribution, with d spatial derivatives, to the induced fermion number is then found by summing over the Matsubara 
poles of the prefactor (p3f), and then using the low temperature expansion in (B21). This gives the following 



(N) 



(d) 

T 



mT 



1 



7T d\{2T) d 



d 2 x 



g(g) 

\v(x)\ 



(d l n a d l n a + \v(x)\)* - {d^d.fi 11 - \v{x)\Y 



(34) 



Note that at any order d of the derivative expansion, the prefactor term in ( p3[ ) leads to a temperature depndent 
prefactor in the induced fermion number that vanishes as T — * 0. Furthermore, the leading low T terms in (|^) have 
a simple form that can be resummed to all orders in the derivative expansion. Combining with the 2-derivative term 
in (EQ), we obtain the following all-orders result for the low temperature correction to the induced fermion number in 
the (2 + 1) dimensional nonlinear a model: 



77} T 1 m 

(N) T ~ (N)o e-9 

7T 



d 2 x 



g(g) 

\m\ 



cosh 



yjdifi a dih a + \v(x)\ 



2T 



cosh 



yjdih a difi a — \v{x)\ 



2T 



(35) 



This is the main result of this paper. It shows explicitly that the topological zero temperature fermion number, given 
by the winding number {N)q in (||) acquires a T dependent correction which is nontopological. This temperature 
dependent correction cannot be expressed in terms of the winding number, and is sensitive to the detailed shape of the 
background field h(x), not just its asymptotic values. Also note that the T — > limit is well defined in the derivative 
expansion regime where ( |l2| ) is satisfied. This resummed expression ([35]) is analogous to a similar all-orders resummed 
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result found in |n| for the (1 + 1) dimensional sigma model background. We find it remarkable that this all-orders 
resummation can also be done for the (2 + 1) dimensional sigma model, as it is a significantly more complicated model. 

Thus far, the only thing we have assumed about n(x) is that it takes values in S 2 , see (0), and that it has a "gentle" 
spatial profile, see (|12|). We now consider two specific explicit forms for the background h(x). 

First, suppose n(x) has the form of a 2-dimensional CV 1 instanton, in which case the winding number (iV)o in (j^) 
is just the instanton number, as is well known for the zero temperature system The instanton field fi(x) satisfies 
the first order instanton equation 



d l h a = ±eije ahc h b djh c 

where the ± signs refer to the anti-instanton/instanton cases respectively (one can easily check this by substituting (|3 

J 3 



(36) 



in the expression (M) to find the sign of (N) ). Now, in general we have that |w(x)| 2 = 2(difi a difi a ) 2 — 2{difi - djh) 2 . But 



for an instanton background n satisfying (|36|), it follows that (dih-djh) 2 = j{difi-difi) 2 . Therefore, \v(x)\ = (dih-difi), 



(37) 



and q(x) = +(9^n • difi). Then the low T expansion ( p5| ) simplifies somewhat 

y / 2dih ■ difi 



, , ITlT 

(N) ± e-T 

7T 



d z x 



cosh 



2T 



To be even more explicit, we choose h to be a fc-instanton solution, in which case (N)q = k. The explicit CV 1 
instantons are most easily expressed in terms of the following stereographically projected fields w 1 and ui . 



2U 1 



2n 2 



(1-n 3 ) ' 

Defining uj = uj 1 + iuj 2 , the original sigma model field h{x) is 



(1-n 3 ) 



(38) 



h(x) 



Keui 
Imcj 



1 



(39) 



1 



As is well known, the instanton equation 



then becomes a simple Cauchy-Riemann condition for lu, so that 



instantons are characterized by a meromorphic function of z 

~(Z - Zq) 



uj(z) = 



x + ix . For example, 

k 



A 



represents k instantons at the location zq, each with a common length scale A. For such an instanton, 

1 



difi ■ difi = 



2fc2 
A 5 * 7 



z - z \ 



\2k-2 



2k\ 



(40) 



(41) 



The spatial integrals in (|37j) can now be done, and one finds 

1 / 2k 



(N) T ~k±—e- 



E 



^ (2ti + 2)! \2i\T 



n=0 



2n 



k' 



(42) 



where f3(m, n) — ^l^r^n) ^ s ^ ne Euler beta function. The first term, k, in ( ^2| ) is the zero temperature fermion 
number (N)o, which is simply the instanton number fc, and which is manifestly independent of the instanton scale 
A. On the other hand, the finite T correction in ( f42| ) is manifestly dependent on the instanton scale A, reflecting its 
nontopological nature. It is interesting to note that this correction term is a convergent series in for any instanton 
number k. For a single instanton (k = 1), it simplifies further to 



A 00 



1 - (4mA)e 



(2n + l)(2n + 2)! \XT 



™. f XT cosh y - 1 
V 2 



dy 



(43) 
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The second physically interesting form for the sigma model background held n is the "baby Skyrmion" case 

cos</>sin#(r) 

n(r,<f>)= \ sm<f>sm6(r) j (44) 
cos 9(r) 

Here, r and <j) are the 2-dimensional polar coordinates and 9(r) is some radial prohle function. The "baby Skyrmion" 
is the 2 dimensional analog of the 3 dimensional Skyrme model for baryons, with a "hedgehog" ansatz for the meson 
fields jl6[ ]. For such an ansatz, the winding number density is 

q(x) = e abc e ij n a d t n b &n c = - sin 9{r) 9' (r) (45) 

r 

so that the zero temperature induced fermion number (^) is 

1 



roc jn 

4tt / dr — sm9(r) 
Jo dr 



L cos(9(0)-cos6l(oo)] (46) 



This shows clearly that (N)o is topological : it only depends on the asymptotic values of the profile function 9(r), not 
on its specific spatial profile or scale. The low temperature correction to the fermion number in this case is obtained 
from (|3g) by noting that for the ansatz 



2 



sin 9{r) x 2 



din a din a = (9'(r)) 
(dS-djh) 2 = (^(r))* + f ^2 V (47) 



in which case, 



|%)0)| = W^M ( 48 ) 
r 



Thus, the low temperature induced fermion number for the "baby Skyrmion" background is 



(N) T ~ <iV) - 2mTe-^ / r dr sign (9'(r) sin 9(r)) 



(\6>(r)\ + ^(111 \ 
cosh — — ■ — cosh 



2T \ 2T 



(49) 



The background field 9{r) has to satisfy certain boundary conditions. Requiring that, for a localized soliton solution, 
the chiral field as r — > oo should equal its vacuum value, i.e, h(r, 4>) should be independent of the angle <fi, fixes 
9{oo) = kit, where k is an integer. At the origin, finitencss of the soliton energy requires 9(0) — nir, where n is an 
integer. As an example, consider the following ansatz for the radial profile function 9(r): 



Then 



Substituting this in (f49|) we obtain 



9(r) = 2 arctan ( ^- ) (50) 
2 A 



/.V), - I - \n,\r - I " dy C ° Shl i 1 (52) 



which is, incidentally, exactly the finite temperature fermion number ( |43| ) for the CV 1 single instanton background. 

The most important thing to be learned from the baby Skyrmion result ( f49| ) is that while the zero temperature 
fermion number (N)q only depends on the asymptotic values of the radial profile function 9(r), the finite T correction 
is sensitive to the details of the shape of 9(r). This means, for example, that in a variational calculation where the 
asymptotic values of 9(r) are fixed (to keep the Skyrme number fixed), but the shape of 9(r) is varied, the finite 
temperature fermion number will vary. 
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IV. CONCLUSIONS 



To conclude, we have computed the finite temperature induced fermion number for fermions in a static nonlinear 
cr-model background in (2 + 1) dimensions. This calculation illustrates the splitting of the induced fermion number 
into a zero temperature piece, which is topological (here, the winding number of the sigma model background), and a 
finite temperature correction which is nontopological p0| . The calculation was done in the derivative expansion limit 
where the spatial derivatives of the background fields are assumed to be small compared to the fermion mass scale 
parameter m. We found that it is possible to resum the derivative expansion to all orders, in the low temperature 
limit where T -C m. Such a resummation was done previously Jl9| in the (1 + 1) dimensional cr-model case, but we 
find it remarkable that it can also be done in the more complicated (2 + 1) dimensional case. These results were then 
applied to two specific background fields in the cr-model case : a 2-dimensional CV 1 instanton, and a 2-dimensional 
"baby Skyrmion" . In each case, the finite temperature correction can be computed in the low temperature limit. 
For an instanton background the final expression for the induced fermion number (N)t becomes a simple convergent 
series in where A is the length scale of the instanton, thus explicitly illustrating the nontopological character of 
the finite temperature corrections in (N)t- 

Finally, we make several comments on possible extensions of this work. As has been emphasized recently at zero 
temperature in |3l], the induced charge (and, more generally, induced current) is only one part of the induced 
effective action for fermions in a cr-modcl background. There is, in addition, a Hopf term which is topological, but 
which does not contribute to the induced current (since this is defined through a functional variation of the effective 
action with respect to a gauge field). This Hopf term has important implications for the spin and statistics of solitons 
|32] , [34| . It would be interesting to investigate these issues at finite temperature, in the language of finite T effective 
actions. Another result from our calculation is the implication it has for the chiral sigma models in (3 + 1) dimensions. 
There, the zero temperature calculation is very similar to the (2 + 1) dimensional case. At finite temperature we learn 
from our calculation that there will be a finite temperature correction and that it will be nontopological - that is, 
while the zero temperature induced fermion number is the Skyrme number, the temperature dependent correction will 
be a much more complicated functional of the Skyrme background. This is already clear from the general argument 
about finite T fermion number in pp| , together with the numerical results (see, e.g., [ 
of the fermion spectrum on the scale of the Skyrme background. That is, since the single-particle fermion spectrum 
is not symmetric and is sensitive to the scale of the background, the finite temperature correction is also sensitive to 
this scale, and thus must be nontopological. From a calculational standpoint, we learn that while at zero temperature 
the induced fermion number arises from the leading nontrivial order of the derivative expansion jl5|, [36j , at finite 
temperature it is necessary to go beyond this order. In other words, at finite temperature there are contributions 
that are not total derivatives, and these contribute to the integrated charge. It is a much more difficult problem to 
investigate higher orders of the derivative expansion in the (3 + 1) dimensional Skyrme model. However, the results of 
this paper for the (2 + 1) dimensional model suggest that it should be possible to explore higher order of the derivative 
expansion, at least in the low temperature limit. It would also be of interest in this context to extend the systematic 
renormalization approach developed in jlj, ^], which is ideally suited to numerical evaluation of induced charges, to 
finite temperature. 



38 ) for the sensitive dependence 
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